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We investigate magnetic Taylor-Couette flow in the presence of an imposed axial mag- 
netic field. First we calculate nonlinear steady axisymmetric solutions and determine 
how their strength depends on the applied magnetic field. Then we perturb these so- 
lutions to find the critical Reynolds numbers for the appearance of wavy modes, and 
the related wavespeeds, at increasing magnetic field strength. We find that values of im- 
posed magnetic field which alter only slightly the transition from circular-Couette flow 
to Taylor-vortex flow, can shift the transition from Taylor-vortex flow to wavy modes 
' by a substantial amount. The results arc compared against onset in the absence of a 

. magnetic field. 
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1. Introduction 



Taylor-Couette flow is one of the most important examples of fluid systems which 
exhibits the spontaneous formation of increasingly complex dynamic flow structures. 
This occurs through a sequence of transitions which take place as the drive is increased 
(Andereck et al. 1986). The first transition takes place when circular-Couette flow (CCF) 
becomes unstable to axisymmetric disturbances which grow into a toroidal flow pattern 
(Taylor- vortex flow — TVF). The second transition occurs when the Taylor vortices 
become unstable to non-axisymmetric perturbations, which result in a variety of time- 



dependent flows (wavy modes). More transitions take place at higher Reynolds numbers 



and a great number of variations of the original pattern observed by Taylor (1923) have 
been investigated (Egbers & Pfister 2000). 

Our concern is Taylor-Couette flow in the presence of an axially imposed magnetic 
field. Despite early theoretical (Chandrasekhar 1961) and experimental interest (Donnelly 
& Ozima 1962), this case has been much less studied. The success of recent dynamo 
experiments (Gailitis et al. 2001; Steiglitz & Miiller 2001) and the astrophysical interest 
in the magneto-rotational instability (Ji et al. 2001; Riidiger & Zhang 2001) add further 
motivations. 

In our previous paper (Willis & Barenghi 2002) we presented a convenient numerical 
formulation of the hydromagnetic Taylor-Couette flow problem. We determined critical 
Reynolds numbers for the appearance of axisymmetric Taylor vortices, calculated finite 
amplitude solutions and compared our results with existing theoretical work and experi- 
ments. The aim of this paper is to continue this investigation into the three-dimensional 
time-dependent flow regime of the wavy modes. 
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2. Formulation and method of solution 

Throughout this work we shall use the same notation as that of our previous paper. 
The fluid is contained between two concentric cylinders of inner radius R\ and outer 
radius i?2 which rotate at constant angular velocity and fl 2 respectively. We make 
the usual assumption that the cylinders have infinite height, use cylindrical coordinates 
(r,8,z), and assume that a constant magnetic field Bo = poHz is applied externally in 
the axial direction. 

We introduce the following dimensionless parameters: radius ratio (77), Reynolds num- 
bers (Rei and Re 2 ), Hartmann number (Q) and magnetic Prandtl number (£) defined 

as 

V = Rl /R 2 , i^^M, i = l,2, Q = ^^! i= V (2.1) 
v pv A 

where S is the gap width, p the density, v the kinematic viscosity, A the magnetic dif- 

fusivity and po the magnetic permeability. Hereafter we assume that p, v, A and po are 

constant. 

The dimensionless equations governing incompressible hydromagnetic flow are 

d t u+{u- V)u = -Vp+V 2 «+|(VAB) AB, V • u = 0, (2.2a, b) 

9 t B = |v 2 B + VA(MAB), V • B = 0. (2.2c, d) 

These equations have as a steady-state solution, it, the circular-Couette flow, where 

u r = u z = 0, u e = Ar + B/r. (2.3) 

The constants A and B are determined by the no-slip boundary conditions. 

Boundary conditions for the magnetic field can have a huge influence on the flow, as 
seen by Hollerbach fc Skinner (2001) when studying the analogous problem in spherical 
geometry. Here the axially imposed magnetic field does not penetrate the boundaries 
but there are still differences between the theoretical results of |Chandrasckhar (1961)| 
for insulating and perfectly conducting cylinders. However, in experiments by Donnelly 
fc Ozima (1962) using mercury, only a small difference was found between results with 
Perspex and stainless-steel containers. Therefore for practical reasons we take insulating 
boundary conditions for the magnetic field. If the magnetic field is expanded over modes 
of the form 

B(r, 9, z) = B(r) S az+me \ (2.4) 
then the boundary conditions are 



a = m = : B e =B z =0: 
a = 0,m/0: dgB r = ±mBo, B z — 0; 

a^0: d z B r = ^^B z , - d e B z = d z B e 



(2.5) 



where for ± we take + at R%, — at i?2- The function B m (r) denotes the modified Bessel 
functions I m (ar) and K m (ar), and the boundary condition ( |2.5| ) is evaluated at Ri and 
i? 2 respectively. 

The magnetic Prandtl number, £, is very small in liquid metals available in the lab- 
oratory (mercury, f = 0.145 x 10~ 6 ; liquid sodium at 120°C, £ = 0.89 x 10~ 5 ). We 
set 

B = B +&, (2.6) 
where Bq is an externally applied field. In the limit £ — > the equations governing the 
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disturbance to the circular-Couette flow become 

{d t - V 2 )w' = N - Vp', V • u' = 0, (2.7a, 6) 

V 2 b = N B , V-6 = 0, (2.7c, d) 

where 

iV = Q(V A 6) A B - (u ■ V)u' - (u' ■ V)u, N B = -V A (u A B ). (2.7e,/) 

The finite amplitude disturbance u' = u — u satisfies homogeneous Dirichlet boundary 

conditions. The magnetic field b satisfies the same boundary conditions as B and is 
completely defined by the velocity at any particular time. 



The method of solution of (2.7 a-d) was the topic of our earlier paper, Willis fc Barenghi 
(2002) , to which we refer for further details, including tests of the numerical method and 
comparison with calculations performed by others. Here it suffices to say that to ensure 
divergence-free fields we adopt the toroidal-poloidal decomposition 

A = %l) 6 + (t)oz + W A (V>r) + V A V A {<j>r), (2.8) 

where ip(r,t,z), <p(r,t,z) and tpo( r ), <fio{ r ) contain the periodic and non-periodic parts 
of the field respectively. The potentials tjj, <j) are expanded over Fourier modes in the 
periodic coordinates and Chebyshev polynomials in the radial direction. 

The governing equations for the magnetic field are the r-components of the induction 
equation and its first curl. Although it is commonplace to take the first and second curls 
for the velocity, we follow the procedure applied to the magnetic field and take the r- 
components of the momentum equation and its first curl. As the pressure has not been 
eliminated, we also take the divergence to obtain the pressure-Poisson equation. 

Whilst the decomposition ( |2.8| ) raises the order of the equations, the extra derivatives 
appear in the periodic coordinates, so no extra boundary conditions are required. All five 
governing equations are second order in r. This property makes them easy to timestep 
stably and accurately, even for fully nonlinear three-dimensional flows, and simplifies 
implementation a great deal. 



The equations are collocated in the radial direction and the nonlinear terms (2.7 ' e,f) are 
evaluated pseudo-spectrally where necessary. Explicit Adams-Bashforth timestepping is 
used on the nonlinear terms and implicit Crank-Nicolson on the linear terms. 



3. An imposed axial field 

The equation defining the magnetic field, (2.7c), may be re-written as 

VA VA6 = VA(tt'Ai), (3.1) 



since V A (u A z) = 0. Equation (3.1) can be integrated immediately to obtain 

VAb=(«'Ai)-V$, (3.2) 
where, by taking the divergence, ^ solves 



V 2 * = t4, for i?i < r < R 2 , 
a r * = 0, onr = i?i,i? 2 , 



(3.3) 



and oj' z is the z-component of vorticity. The function VE" is related to the z-component 
of the streamfunction, for three-dimensional incompressible flow (u = V A 3>; u> = 
— V 2 <fr). The boundary conditions for <1> are determined by the no-slip condition. Here, 4" 
satisfies the homogeneous Neumann boundary condition that there be no current through 
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the boundaries, f ■ V A b — 0. The following analysis may differ slightly for non-insulating 
boundaries. 

The effect of the magnetic field on the fluid is introduced via the Lorentz force, which 
can now be written as 

L = Q(VA6)Az = -Q K, u' e , f - Q [r^dgV, -d r V, of . (3.4) 

The first term on the right hand side is directly proportional and opposed to any flow 
across the imposed magnetic field lines. Axial flow is unaffected by the Lorentz force. 
Therefore it is not surprising that the flow pattern is found to elongate in this direction 
as the strength of the imposed field increases (Chandrasekhar 1961). Taking the curl of 

& 

V A L = Q (d z u A z - oj' z z) - Q (V<9 Z # -u' z z). (3.5) 



The w^-terms on the right hand side of (3.4) cancel. It is well known that components of 



vorticity perpendicular to an imposed field are preferentially damped. The vorticity of the 
underlying circular-Couette flow (CCF) points in the z direction, so CCF is unaffected 
by the imposed magnetic field, and the effective viscosity of the fluid is unchanged by its 
presence until the appearance of the Taylor instability. 



Taking the dot-product of u' with (2.7a) and integrating over the volume gives an 



energy balance equation for the disturbance. The nonlinear term 

J u' ■ (u 1 ■ V)u' dV = | J V ■ (u' 2 u') dV (3.6) 

vanishes due to the boundary conditions. This term represents the transfer of energy 
within the disturbance to higher modes (smaller length scales) by advection. Using the 
property 

J u' • V V dV = J V • (u' A w') dV-Ju'-u' dV, (3.7) 

where similarly the divergence integral vanishes, the change in kinetic energy can be 
written as 

| dt J U' ■ U' dV = J u' r u' g (J- - d r ^j ilg dV 

- J u/ • u>' dV - Q J « 2 + u'g 2 - u' g d r ^) dV. (3.8) 

As the instability that first appears is initially axisymmetric, azimuthal derivatives have 
been ignored in ( |3.8| ). For a steady-state solution the terms on the right hand side must 
balance, hence these three terms can be identified as 

E 1 -E 2 -E z = 0. (3.9) 

The first term, E\, is the supply of energy to the disturbance from radial shear of the 
circular-Couette flow. The second term, E 2 , is the viscous dissipation and the third, E3, 
is the magnetic damping term. 



4. Stability of CCF in a weak magnetic field 

The linear stability of this flow at various radius ratios (Soundalgekar, Ali & Takhar 
1994), for co- and counter-rotating cylinders (Chen & Chang 1998), and at asymptotically 
large Q (Chandrasekhar 1961), is a well studied topic, so we focus briefly on the behaviour 
at small Q. In this and the following sections we find the magnetic field does not need to 



Hydromagnetic Taylor-Couette flow 



5 



be large, Q = 0(10), in order to have a significant affect. By experimental standards this 
Q is not very large; experiments with Q > 10 3 were performed by Donnelly & Ozima 



(1962) 



The disturbance to the circular-Couette flow can be decomposed into Fourier modes, 

oo 

u '( r ? z ) = 75 u o( r ) +XXfc( r ) cosakz ' (4- 1 ) 



fc=i 



where 2ir/a is the critical wavelength at which the fundamental disturbance (k = 1) first 
appears. Due to orthogonality of the cosine function we have 

i [u'.u'dV= -J2 [<-<rdr, (4.2) 

J a k J 

which allows us to consider the contribution of each mode to the energy. 

For a disturbance to the fundamental mode of infinitesimal amplitude A\ , the nonlinear 
term ( |3.6| ) will be very small and modes k > 1 can be assumed negligible. At the critical 
Reynolds number this disturbance neither grows nor decays and can be considered in 
isolation. 

With a fixed outer cylinder the energy, E\, supplied to the fundamental disturbance 
is proportional to Re\, the driving imposed by the inner cylinder. The inviscid Rayleigh 
criterion predicts that the flow is unstable for any rotation of the inner cylinder. The 
viscous dissipation, E2, prevents the immediate onset of TVF. It originates from shears 



within, and so derivatives of, the disturbance. From ( |3.8|) and (4.1) for the first mode, 
it is expected that E2 — 0(a 2 Af). For a weak magnetic field the magnetic damping, 
E3, is 0{QA\), linearly proportional to Q. Provided that the magnetic field is not large 
enough to affect the wavenumber of the flow, we expect that the critical driving, 

R&c j at 

which the instability first appears, to be delayed linearly with Q. Further, since a ss 3 
for hydrodynamic flows, if Q — 0(10) then E2 and E3 are expected to be of comparable 
size. 

Results of calculations at radius ratios 77 = 0.65, 0.72, 0.83 are given in figure |. The 
wavenumber decreases approximately linearly by 15% over the not-so-small range of Q 
in the figure. For each 77 the wavenumber a = 3.13, 2.93, 2.70 ± 0.01 at Q = 0,15,30 
respectively. The delay of transition to TVF remains almost linear in Q. In the limit 
Q^Owe find that Re c (Q) / Re c (0) = 1 + 7Q where 7 w 0.0181 for all 77 in figure Q. As 
no energy is passed to higher modes for very small A\ , the structure of the disturbance 
does not change appreciably with Q or r\. Therefore the ratio of magnetic to viscous 
dissipation remains the same with different r\. 



5. Nonlinear axisymmetric TVF 

In this section we study the nonlinear flows that occur above the first transition. If 
the inner cylinder is driven past the critical rate, then energy is quickly transferred to 
higher modes by the nonlinear advection. Isolating one of these modes, the dissipative 
terms are E 2 = 0{a 2 k 2 A\) and E 3 = 0(7 QA 2 .). Viscous dissipation quickly becomes 
dominant as energy moves to higher modes. This can be seen in figure |^. However, 
magnetic damping remains significant on the fundamental mode, draining energy before it 
can be passed to higher wavenumbers. There the fluid viscosity determines the amplitude 
of the disturbance. The net result can be seen in figure || where we plot the amplitude of 
the disturbance at different values of Q. A relatively wide gap, r\ = 0.65, was chosen for 
the calculations, as the axisymmetric flow is stable to non-axisymmetric perturbations 
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Figure 1 . Delay of onset of TVF by an imposed axial field. Dotted lines are linear 

extrapolations from Q small. 
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Figure 2. Typical spectrum of energy terms. For the first few modes magnetic (dashed) and 
viscous (solid) dissipations are comparable. Viscous dissipation dominates the behaviour of 
higher modes, and E 2 /E 3 oc k 2 . (Q = 30, r) = 0.65, a = 2.71, Rei = 1.5Re c , Re c = 113.4.) 



well beyond the onset of TVF. The maximum axial velocity, rather than the radial 
velocity, was used to characterise the amplitude of nonlinear TVF, as the outflow becomes 
increasingly jet-like as Re\ is increased. Tabeling (1981)| used an amplitude expansion 
about the point of criticality in the narrow gap limit to show that in the weakly nonlinear 
regime A oc P{Re\ — Re c )^ , and that f3 does not strongly depend on Q. From figure || 
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Figure 3. Square of maximum axial velocity versus Rei in the presence of imposed axial 

magnetic fields. r\ — 0.65 . 

we see that this approximation is valid well beyond the critical point, even though the 
nonlinearity affects the flow pattern. 

Kikura, Takeda & Durst (iyy9) performed experiments to measure fluid velocities 
in nonlinear hydromagnetic flow. Their imposed field was different from ours, so we 
cannot make a direct quantitative comparison. However, the same order-of-magnitude 
arguments on the energy of the harmonics still hold, and the numerical results presented 
in this section are qualitatively the same as their experimental results. 

Since the magnetic damping is effective only on larger length scales, apart from a 
change in wavelength of the flow, the flow pattern of nonlinear axisymmetric TVF in 
the presence of an imposed axial field is very similar to that in the hydrodynamic (non- 
magnetic) case. This is apparent from figure |§| where the main difference is only that the 
onset of TVF is delayed. This can also be seen in figure f| where (Rei — Re c (Qj) is the 
same for both plots (with r] = 0.83, in (a) Q — 0, Rei — 1.5Re c (0), a — 3.13 and (6) 
Q = 20, Rei = Re c (Q) + 0.5i?e c (0), a = 2.85). Although the axial wavelengths differ, 
both are plotted the same size for comparison. The azimuthal flow is faster than circular- 
Couette flow (u' g > 0) at the outflow regions (z — ir/a) and slower (u' e < 0) at the inflow 
regions (z — 0, 2ir/a). The only clear observable difference is that the fast azimuthal jet- 
flow region that occurs at the outflow is slightly larger for the magnetic flow, in addition 
to the stretched axial wavelength. However, figure ^& is almost indistinguishable from 
hydrodynamic flow at the longer wavelength, a = 2.85, where Rei = 1.5Re c (not plotted). 
The deviation from circular-Couette flow is large; the maximum of u' g is approximately 
35% of the maximum of ug — Rei. The maximum values of lu^ezl for all three sets 
of parameters are within 2%, except for \u' z \ in the magnetic flow which is 8% larger. 
The Lorentz force has no z-component and so does not damp axial flow, and therefore 
the curves in figure || are not quite parallel. Despite only small differences in the flow 
patterns, only the magnetic flow is stable to non-axisymmetric disturbances, which is the 
topic of the next section. 
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Q = 20, Rex = Re c (Q) + 0.5Re c (0), a = 2.85. 



6. The transition to wavy flow 

Finding the preferred mode and transitions between fully developed nonlinear wavy- 
modes is an enormous challenge. For the sake of simplicity, we investigate which modes 
are possible by monitoring the growth or decay of infinitesimal non-axisymmetric distur- 
bances to the axisymmetric TVF. The disturbance translates in the azimuthal direction 
at some fraction of the rotation rate of the inner cylinder — the wavespeed s. 

In the hydrodynamic case, Jones (1985) calculated stability boundaries and corre- 
sponding wavespeeds by solving the eigenvalue problem. The stability of axisymmetric 
flow was found to depend strongly on the radius ratio. Figure]^ shows the stability bound- 
aries and corresponding wavespeeds at onset of the wavy modes. The stability boundaries 
are shown relative to Re c for the onset of TVF (m = 0). For increasing Reynolds num- 
bers and n > 0.75, we see that TVF can be destabilised to the m = 1 wavy mode (lower 
part of the boundary in figure ||a) and then be restabilised by a further increase of the 
Reynolds number (upper part of the boundar y). Such transition seq uences predicted by 



Jones' results were verified experimentally by Park fc Jeong (1984) 



The critical Reynolds number for the disappearance of m = 1,2 modes (upper bound- 
ary) increases rather rapidly as r\ is increased past about 0.85. Therefore the choice 
j] = 0.83 was taken for our magnetic calculations. 

The magnetic field, imposed axially, has the effect of increasing the critical wavelength 
for onset of TVF. Jones (1985) held a fixed at 3.13, but this becomes inappropriate for our 
magnetic calculations. Figure | shows the stability of the equilibrated TVF that appears 
at the critical wavenumber for the given field strength, a(Q). The dotted lines are the 



stability boundaries for bifurcation from circular-Couette flow, as calculated by Chen & 



Chang (1998). The difference between the dotted lines and the solid lines confirm Jones' 



assertion that it is essential to perturb TVF (which must be computed numerically) to 
understand the onset of wavy modes. The stability boundaries for the onset of wavy 
modes are sensitive to small differences in the parameters, and the magnetic field does 
not need to be strong to produce interesting effects. 
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Figure 5. (a) Stability of axisymmetric TVF to non-axisymmetric perturbations over a range 
of 7], as determined by Jones (1985), a — 3.13. (b) Corresponding wavespeeds, s, at onset as a 
fraction of fii. The dotted line is the intersection with the magnetic results at rj — 0.83 . 



The magnetic damping term in ( p.8| ) damps the radial and azimuthal components of 
the axisymmetric flow. The curvature of the cylinders also has an additional effect on the 
viscous dissipation of these two components of the velocity. There is a clear similarity 
between increasing the imposed magnetic field strength and decreasing the radius ratio. 
We do not pursue this relationship too far as the stability boundaries have not been fully 
explained, even for hydrodynamic flows. However, significant progress on the mechanism 



involved was made by Jones (1985) 



Jones (1985) noted that dissipation due to the radial and axial shear within TVF 
appear to be important factors in the stability of the azimuthal jet that appears at 
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Figure 6. For r\ = 0.83, a = a c {Q) (a) stability of hydromagnetic TVF to non-axisymmetric 
perturbations for increasing magnetic field strength, and (b) corresponding wavespeeds at onset. 
Dotted lines are bifurcations from CCF. 



the in- and out-flows. The elongation of Taylor cells over the range of Q in figure [3 
(by around 10% ) might be a factor in the suppression of wavy modes. [Antonijoan fc 
Sanchez (2002) investigated the hydrodynamic stability to wavy modes as a function of 
a. The resulting stability boundaries are qualitatively similar to figure |6|a for increasing 
wavelength. However, from the results of Antonijoan fc Sanchez (2002)| , the change in a 
over the range for Q in figure is a factor of 2-3 smaller than that required to suppress 
only the m = 1 mode. The enhanced stability must also be attributed to magnetic 
damping of the disturbance. 

The general trend in figure ||& is for the wavespeed, s, to decrease with decreasing r\. 
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For wider gaps there is more fluid in the slower outer regions and the mean azimuthal 
flow speed decreases with increasing gap width, as does s. Also, as the Reynolds number 
is increased the outflow region narrows and the inflow region becomes larger than that 
of the outflow. Therefore the bulk of the fluid is in the slower inflow region and again 
the mean azimuthal flow is reduced. This mean-flow interpretation of the wavespeed is 
consistent with the findings in figure |6]&. Here the wavespeed decreases with increasing 
Reynolds number, as the outflow still narrows and the inflow broadens, but the gap width 
is fixed and there is no additional drop in the pattern when the magnetic field strength 
increases. Compare figure fyb for decreasing ij with figure ^& for increasing Q. 



7. Conclusions 

We have investigated magnetic Taylor-Couette flow in the nonlinear regime. In the 
presence of an imposed axial magnetic field the Lorentz force is found to have a significant 
damping affect, but only at larger length scales. This is a consequence of the small 
magnetic Prandtl number limit, relevant to experiments with liquid metals, for which 
the magnetic field is completely defined by the velocity field. 

When the imposed field is not too large, the axial wavenumber of the flow is not greatly 
affected. In this weak-field regime the stability of circular-Couette flow is enhanced lin- 
early with Q, Re c (Q) / Re c (Q) = 1 + -fQ where 7 appears to be independent of the gap 
width. 

In the nonlinear regime, magnetic damping affects disturbances at the fundamental 
axial wavenumber, but the remaining energy passed to higher modes is dissipated mainly 
by the fluid viscosity. This determines the amplitude of the disturbance, which behaves 
approximately like j3(Re\ — Re c (Q))^ well beyond the critical point, where j3 does not 
depend strongly on Q. This is consistent with the findings of the amplitude expansion 



by Tabeling (1981 



Our main finding is that the magnetic field has its most striking effect on the stability 
of TVF to wavy modes. A small field is capable of pushing the secondary instability from 
only a few percent above the first instability to several times past the critical Reynolds 
number for the onset of TVF. This is similar to the relative stability of TVF to wavy 
perturbations in wide gaps. As in hydrodynamic flows, the wavespeed decreases with 
increased Reynolds number, but the dependence of the wavespeed on the imposed field 
strength does not appear as strong as the dependence on the gap width. 

The significant enhanced stability observed in the calculations above occurs at only 
relatively small imposed field strengths, well within experimental range. As a secondary 
bifurcation, the transition to wavy modes is difficult to detect accurately via torque 
measurements. A visualisation technique, using ultrasound, is being developed for opaque 



fluids by Kikura, Takeda & Durst (1999). This type of transition could serve as a good 
test for flow visualisation in magnetic fluids, where the transitions are accurately defined. 
The new technique is of particular interest to those working on dynamo experiments. 

The authors wish to thank Prof. Chris Jones for helpful comments and suggestions 
during this work. 
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